Experimental evidence of a phase transition in a closed turbulent flow 
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We experimentally study the susceptibility to symmetry breaking of a closed turbulent von 
Karman swirling flow from Re = 150 to Re ~ 10 6 . We report a divergence of this susceptibil- 
ity at an intermediate Reynolds number Re = Re x ~ 90 000 which gives experimental evidence that 
such a highly space and time fluctuating system can undergo a "phase transition" . This transition 
is furthermore associated with a peak in the amplitude of fluctuations of the instantaneous flow 
symmetry corresponding to intermittencies between spontaneously symmetry breaking metastable 
states. 
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Phase transitions are ubiquitous in physical systems 
and generally are associated to symmetry breakings. For 
example, ferromagnetic systems are well known to un- 
dergo a phase transition from paramagnetism to ferro- 
magnetism at the Curie temperature T c . This tran- 
sition is associated with a symmetry breaking from 
the disordered paramagnetic — associated to a zero 
magnetization — toward the ordered ferromagnetic phase 
— associated to a finite magnetization — [1] . In the vicin- 
ity of T c , a singular behaviour characterized by criti- 
cal exponents is observed, e.g. for the magnetic sus- 
ceptibility to an external magnetic field. In the con- 
text of fluid dynamics, symmetry breaking also governs 
the transition to turbulence, that usually proceeds, as 
the Reynolds number Re increases, through a sequence 
of bifurcations breaking successively the various symme- 
tries allowed by the Navier-Stokes equations coupled to 
the boundary conditions [2 . Finally, at large Reynolds 
number, when the fully developed turbulent regime is 
reached, it is commonly admitted that all the broken 
symmetries are restored in a statistical sense, the statis- 
tical properties of the flow not depending anymore on Re 
[3]. However, recent experimental studies of turbulent 
flows have disturbed this vision raising intriguing fea- 
tures such as finite lifetime turbulence ^ — questioning 
the stability of the turbulent regime — and possible exis- 
tence of turbulent transitions [SHTTj. Consequently, de- 
spite turbulent flows are intrinsically out-of-equilibrium 
systems, one may wonder whether the observed transi- 
tions can be interpreted in terms of phase transitions 
with a symmetry-breaking or susceptibility divergence 
signature. In this paper, we introduce a susceptibil- 
ity to symmetry breaking in a von Karman turbulent 
flow and investigate its evolution as Re increases from 
150 to 10 6 using stereoscopic Particle Image Velocime- 
try (PIV). We observe a divergence of susceptibility at 
a critical Reynolds number Re — Re x ~ 90 000 which 
sets the threshold for a possible turbulent "phase tran- 
sition". Moreover, this divergence is associated with a 
peak in the amplitude of the fluctuations of the flow in- 
stantaneous symmetry. 




FIG. 1. Schematic view of the experimental setup and the 
impellers blade profile. The arrow on the shaft indicates the 
impeller rotation sense studied. Symmetry : the system is 
symmetric regarding any 7?-n--rotation of angle n around any 
line of the equatorial plane which crosses the rotation axis. 

Our experimental setup consists of a Plexiglas cylin- 
der of radius R = 100 mm filled up with either water or 
water-glycerol mixtures. The fluid is mechanically stirred 
by a pair of coaxial impellers rotating in opposite sense 
(Fig. [lj. The impellers are flat disks of radius 0.925 R, 
fitted with 16 radial blades of height 0.2 R and curvature 
radius 0.4625 R. The disks inner surfaces are 1.8 i? apart 
setting the axial distance between impellers from blades 
to blades to 1.4 R. The impellers rotate, with the con- 
vex face of the blades pushing the fluid forward, driven 
by two independent brushless 1.8 kW motors. The rota- 
tion frequencies j\ and f% can be varied independently 
from 1 to 12 Hz. Velocity measurements are performed 
with a stereoscopic PIV system provided by DANTEC 
Dynamics. The data provide the radial u r , axial u z and 
azimuthal u v velocity components in a meridian plane 
on a 95x66 points grid with 2.08 mm spatial resolution 
through time series of 400 to 27 000 fields regularly sam- 
pled, at frequencies from 1 to 15 Hz, depending on the 
turbulence intensity and the related need for statistics. 
The control parameters of the studied von Karman flow 
are the Reynolds number Re — 7r(/i + f2)R 2 /v, where v 
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FIG. 2. (a) and (b) Schematic drawings of the flow topology 
and (c) and (d) corresponding experimental maps of mean ve- 
locity field of the turbulent von Karman flow at Re = 800 000 
for (c) 9 = (7 = 0) and (d) 9 = -0.0147 (7 < 0). The 
color maps the azimuthal velocity u v , from blue to red ("jet" 
colormap), whereas the arrows map the (u r ,u z ) field. The 
resolution has been reduced by a factor 2 for better visibil- 
ity. The r o — r symmetry in (c) and (d) reveals that the 
time-averaged mean fields are axisymmetric. 

is the fluid viscosity, which controls the intensity of tur- 
bulence and the rotation number 6 = (A — /a) /(A + A), 
which controls the asymmetry of the forcing conditions. 
The rotation frequencies A, 2 ar e regulated by servo loop 
control and we obtain for 9 a typical absolute precision 
of l%o and time fluctuations of the order of ±2 x 10~ 4 . 
The correlation between these fluctuations and the flow- 
dynamics are negligible. 

When 9 = 0, the experimental system is symmetric 
with respect to any T^-rotation exchanging the two im- 
pellers: the problem conditions are invariant under n- 
rotation around any radial axis passing through the cen- 
ter of the cylinder (Fig. [I]) . The symmetry group for such 
experimental system is 0(2) [12] . When 9^0, the exper- 
imental system is no more TZ n -symmetric, the symmetry 
switching to the £0(2) group of rotations. However, the 
parameter 9, when small but non-zero, can be considered 
as a measure of the distance to the exact 0(2) symmetry: 
the stricto sensu 50(2) system at small 9 can be consid- 
ered as a slightly broken 0(2) system [13l[T4]. Depend- 
ing on the value of 9, the flow can respond by displaying 
different symmetries: (i) the exact 72.^-symmetric flow 
composed of two toric recirculation cells separated by an 
azimuthal shear layer located at z = when 9 = (Figs. 
[2]ja) and (c)); (ii) an asymmetric two-cells flow, the shear 
layer being closer to the slowest impeller (z / 0), when 
9^0 (Figs. [2jb) and (d)); (hi) and finally, a fully non 
symmetric one-cell flow, the whole shear layer being con- 
centrated in between the blades of the slowest impeller, 
when 9 becomes large enough [TO] US] . 

In order to quantify the distance of the flow to the 
T^-symmetry, we use the normalized and space-averaged 
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FIG. 3. Global angular momentum I(t) as a function of 
time for an experiment performed at Re = 43 000 for 8 = 0. 
Green/light gray dots are PIV data sampled at 15 Hz and 
blue/dark gray dots correspond to 1 Hz low-pass filtered data 
If(t). Eye-guiding lines have been drawn at I(t) = ±J* = 
±0.04. 

angular momentum I(Re, 9, t) as order parameter: 

I(t) = U rdrdtpdz ™f] f , 
v Jv tt-R (A + A) 

where V is the volume of the flow [25]. An example of 
time variation of I(t) at 9 = in the turbulent regime 
is provided in Fig. [3] We assume that ergodicity holds, 
meaning that the instantaneous turbulent flow is explor- 
ing in time its energy landscape according to its statis- 
tical probability. In this framework, the time average 
value I of I(t) is equivalent to a statistical mechanics en- 
semble average providing the average is performed over 
a long enough duration in order to correctly sample the 
slowest time-scales. Then, using this ensemble averaged 
order parameter, we define a susceptibility of the flow to 
symmetry breaking xi as: 

dl 



Note that I is proportional to the mean altitude z s of 
the shear layer which is the natural measure of the flow 
symmetry. Contrary to z s , I(t) is defined for any instan- 
taneous velocity field, including turbulent ones. 

In the non-fluctuating laminar case, when 9 = 0,1 = 
due to the symmetry of the flow. In contrast, as 9 drifts 
away from 0, the value of the angular momentum I be- 
comes more and more remote from zero as the asymmetry 
of the flow grows. In such a framework, there is a formal 
analogy between ferromagnetic and turbulent systems. 
For ferromagnetism (resp. turbulence), the order param- 
eter is the magnetization M(T, h) (resp. the angular mo- 
mentum I(Re,9)); the symmetry breaking parameter is 
the external applied field h (resp. the relative driving 
asymmetry 9); the control parameter is the temperature 
T (resp. the Reynolds number Re, or a function of it). 

In the sequel, we first investigate the influence of tur- 
bulence on / and Xi as R e increases from 150 to 10 6 . In 
the laminar flow at Re = 150, the symmetry parameter 
/ = I(t) evolves linearly with 9 (Fig. [4^c)) and the sus- 
ceptibility is xi — 0.240±0.005. Increasing the Reynolds 
number, one expects to reach fully developed turbulence 
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FIG. 4. Reynolds number dependence of: (a) Susceptibility 
to symmetry-breaking \i 01 the von Karman mean flow at 
8 — and (b) standard deviation aj of the global angular 
momentum I(t) at 9 = 0. In (a), the dotted line corresponds 
to the mean field theory approach with a critical Reynolds 
number Re x — 70 000 and the dashed line to Re x = 90 000; 
(c) global angular momentum / and (d) standard deviation 
<tj as a function of 9 for Re = 150 (green □), Re = 65 000 
(red o) and Re = 800 000 (blue A). In (a) and (b), horizontal 
error bars are of the order of the marker size. In (a) vertical 
error bars are computed as the maximum error that can arise 
from the sharp dependence of oi with 9 around 9 — 0. 



around Re = 10 000 [TO]. In such turbulent regimes, ve- 
locity fields of von Karman flows are characterized by 
a high level of intrinsic fluctuations, i.e. fluctuations 
of the same order of magnitude than the mean values 
[15) . Therefore, even when 9 = 0, the TZ^ -symmetry is 
of course broken for the instantaneous flow. However, 
as usually observed for classical turbulence, this symme- 
try is restored for the time-averag cd flow (Fig. gc)), 
which proves that time averages are long enough to cor- 
rectly sample the slowest flow time scales. Then, as in 
the case of the laminar flow, when 9 is varied, we observe 
the breaking of the IZ^ -symmetry of the mean flow (Fig. 
Id)). 

In Fig. gc), we see that, at Re = 800 000, in the 
close vicinity of 9 = 0, 1(9) evolves actually much more 
rapidly with 9 than in the laminar case with the suscep- 
tibility being larger by more than one order of magni- 
tude: xi — 9 ± 1. Therefore, turbulence seems to en- 
hance dramatically the sensitivity of the flow to symme- 
try breaking. Furthermore, for intermediate Re = 65 000, 
the slope of 1(9) around 9 = is even much steeper - 
Xi = 43 ± 1— than for Re = 800 000. In Fig. ga), we 
plot the susceptibility with respect to Re. We see that 
the susceptibility actually grows by more than two orders 
of magnitude — from 0.24 to 46 — between Re = 150 
and Re ~ 90 000, before decreasing by a factor 4 be- 
tween Re ~ 90 000 and Re = 800 000. These results 
suggest a critical behaviour for x/(i?e) near Re = Re x = 



90 000 ± 10 000: a divergence — revealing a continuous 
second order phase transition — or only a maximum — 
revealing either a subcritical bifurcation or a continuous 
transition with finite-size effect. This cannot be exper- 
imentally tested further since the highest measured Xi 
are already of the order of the highest measurable value 
considering the 9 precision of our setup. For higher \9\, 
we observe a crossover — for the slope — in the curve 1(9) 
at \9 r \ = (6 ± 1) x 10" 3 for Re = 800 000 and very close 
to 6 = 0, at \9 r \ = (0.9 ± 0.15) x 10~ 3 , for Re = 65 000 
(Fig. |4|c)). For \9\ > \8 r \, we recover the laminar slow 
evolution of / with 9 up to 9 = ±0.1 where the flow bifur- 
cates to the one-cell topology (not shown). Since 1(9) is 
quite independent of Re for \9\ > \9 r \ at large Re, we can 
extrapolate this linear behaviour to 9 = 0. This extrapo- 
lation describes the ideal behaviour at critical Reynolds 
number Re x if xi diverges: a jump of / between —Iq 
and +Jq where I — 0.05. This can be interpreted as a 
spontaneous "turbulent momentization" Iq at 9 = - 
possibly affected by finite-size effects — by analogy with 
the spontaneous magnetization Mq at zero external field 
for ferromagnetism. It is also similar to the experimental 
results of [To] . 

A signature of this momentization can be seen on the 
instantaneous global angular momentum I(t) for Re near 
the peak of susceptibility and 9 = (e.g. in Fig. [3]). In- 
deed, one observes that I(t) does not remain near zero 
(its mean value) but shows a tendency to lock pref- 
erentially on the plateaus ±7* with 7* = 0.04 ~ To- 
Therefore, the turbulent flow explores a continuum of 
metastable symmetry breaking patterns evidenced by 
-0.04 < I f (t) < 0.04, I f (t) being the 1 Hz low-pass 
filtered value of I(t). The global angular momentum ac- 
tually fluctuates very much along time with two sepa- 
rate time scales: fast fluctuations related to "traditional" 
small scale turbulence and time intermittencies corre- 
sponding to residence time of few tens of seconds. If 
one performs a time average over one of these intermit- 
tent periods only, one obtains a time localized "mean" 
flow, which breaks spontaneously the symmetry, analo- 
gous to what is obtained for true mean flows when 9 =/= 
as presented in Figs. [2jb) and (d). 

The presence of strong fluctuations is not surprising 
here: close to a phase transition we expect critical fluc- 
tuations. To check this, we compute the standard devia- 
tion a i (Re, 9). Fig. gb) shows, for 9 = 0, how 07 varies 
from zero in laminar case to finite values for highly turbu- 
lent flows going through a maximum at Reynolds number 
Re = Re a = 45 000 ± 10 000 located below the peak of 
susceptibility at Re x . Additionally, in Fig. gd), the 
dependence of 07 as a function of the symmetry control 
parameter 9 reveals a strong difference between the two 
Reynolds numbers shown: crj(9) presents a sharp and 
narrow peak at 9 = for Re = 65 000, which does not 
exist for Re = 800 000. The amplitude of this peak from 
its bottom to its top actually measures the additional 
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amount of low frequency symmetry fluctuations due to 
the multistability. This amount of fluctuations appears 
to be connected to the susceptibility increase below Re x . 

The previous experimental results set a strong connec- 
tion between the spontaneous symmetry fluctuations of 
the flow and the mean flow response to the system sym- 
metry breaking: the interpretation of the large fluctua- 
tions of I(t) in terms of multistability suggests that the 
strong observed linear response of the mean flow (Fig. 
gc)) with respect to 9 in the close vicinity of 6 = is 
the result of a temporal mixing between the metastable 
states in different proportions. 

Nevertheless, as the Reynolds number is varied, two 
distinct maxima have been evidenced in our turbulent 
flow: one for the susceptibility xi to the 0(2)-symmetry 
breaking near Re = Re x ~ 90 000, and the other for 
the standard deviation of the global angular momentum 
I(t) near Re = Re a ~ 45 000. In this Reynolds number 
range, the turbulence is generally expected to be fully de- 
veloped, i.e. any non-dimensional characteristic quantity 
of the flow should be Re- independent. This is definitively 
not the case in this von Karman experiment. Actually, 
visual observations of the flow reveal an increase of the 
average azimuthal number m of large scale vortices in the 
shear-layer from m — 3 to m — 4 through an Eckhaus- 
type transition, between Re a and Re x . In the following, 
we make the hypothesis that at least one critical phe- 
nomenon exists in the range 50 000 < Re < 100 000. 
Using the turbulent von Karman-ferromagnetism anal- 
ogy, we can check how classical mean field predictions 
for second order phase transition apply to our system. 
In terms of susceptibilities, it predicts a critical diver- 
gence at a temperature T c : x oc |T — T c | . Since the 
logarithm of the Reynolds number has already been pro- 
posed as the control parameter governing the statistical 
temperature of turbulent flows — T ~ l/logi?e [T7] — , 
this prediction translates in our case into: 

Xi oc 1 1/ log i?e - l/logi?e c | _1 . 

This formula reasonably describes our data with Re c be- 
tween typically 70 000 and 90 000 (Fig. ga)) support- 
ing the asymmetry between the two branches even with 
a unique exponent —1. As far as fluctuations are con- 
cerned, it is difficult to find a reasonable critical exponent 
for 07. However, the results show that the maxima for 
Xi and 07 are clearly separated. This is at variance with 
classical phase transition theory, stating e.g. in the Fluc- 
tuation Dissipation Theorem that aj should be propor- 
tional to xi- A reason for this discrepancy lies in the high 
level of intrinsic fluctuations in our system that makes 
this transition non-classical. Instabilities or bifurcations 
occurring on highly fluctuating systems are commonly 
found in natural systems and the literature reports tran- 
sitions and symmetry breaking at high Reynolds num- 
ber (see, e.g., references [5 HTT1 ITtJI ITM2"U] ) but the corre- 
sponding theoretical tools are still today not well settled. 



Existing studies of phase transitions in the presence of 
fluctuations generally considers systems in which an ex- 
ternal noise — additive or multiplicative — is introduced 
[2 1) . In particular, it has been shown in models that mul- 
tiplicative noise can produce an ordered symmetry break- 
ing state through a non equilibrium phase transition [22] . 
This behaviour could be at the origin of our observed 
transition. Finally, we can notice that, as the Reynolds 
number increases, the "turbulent momentization"/ first 
increases and then decreases, contrary to the magneti- 
zation in the usual para-ferromagnetic transition. This 
result is reminiscent of a reentrant noise-induced phase 
transition similar to that observed in the annealed Ising 
model [23j [24]. The study of the evolution of Iq and/or 
7* with Re requires more statistics and is left for fu- 
ture work. Finally, our turbulent system, in which we 
can have access both to the spatiotemporal evolution of 
the states and to the mean thermodynamic variables, ap- 
pears as a a unique tool to study out-of-equilibrium phase 
transitions in strongly fluctuating systems. 

We thank K. Mallick for fruitful discussions. PPC was 
supported by Triangle de la Physique. 
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